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We study, at zero temperature, the dynamics of a sin- 
gle electron in a Holstein model augmented by site-diagonal, 
binary-alloy type disorder. The average over the phonon 
vacuum and the alloy configurations is performed within a 
generalized dynamical coherent potential approximation. We 
present numerical results for a Bethe lattice with infinite co- 
ordination number. In particular, we investigate, in the in- 
termediate electron-phonon coupling regime, the spectral and 
diffusion properties in the vicinity of the high-energy edge of 
the lowest polaronic subband. To characterize the diffusion 
properties, we define a spectrally resolved delocalization time, 
which is, for a given energy, the characteristic time scale on 
which the electron leaves a given site. We find the delocal- 
ization times substantially enhanced for states with a large 
phonon content, i.e., in the absence (presence) of alloy-type 
disorder at the high-energy edge(s) of the polaronic subband 
(mini-subbands). According to their delocalization times, we 
discriminate between "fast" quasi-particle-like and "sluggish" 
defect-like polaron states and qualitatively address the issue 
of trapping of an electronic carrier. 

PACS numbers: 71.38. -fi, 72.10.Di, 71.35.Aa 



I. INTRODUCTION 

In nature there is an abundance of "polaronic" com- 
pounds in which, due to strong electron-phonon coupling, 
electrons and phonons lose, in a certain temperature and 
density range, their individual identity and new enti- 
ties emerge: polarons jl]. Examples of current inter- 
est are, among others, the high-temperature supercon- 
ducting perovskites 1^, the "colossal" magneto-resistance 
manganites ||], and the titanates The details of 

the actual materials notwithstanding, it is commonly as- 
sumed that, at least qualitatively, many aspects of these 
compounds can be understood on the basis of the Hol- 
stein model supplemented by various terms, such as, 
e.g., Coulomb repulsion, spin-orbit coupling, and specific 
forms of material imperfections. 

In its basic form the Holstein model comprises a single 
tight-binding band locally coupled to dispersionless op- 
tical phonons. Despite this simplicity, the properties of 
the Holstein model are not yet fully understood. Even in 
the extreme dilute limit, that is the case of a single po- 
laron in an otherwise empty band, many basic questions 



remain, despite impressive theoretical progress made in 
recent years p|-p^. 

In particular, the notoriously difficult issue of trap- 
ping pO|-p^ is far from being completely resolved. The 
trapping properties are only understood for the polaron 
groundstate of the ordered Holstein model. In an im- 
portant paper, Lowen [Tc| ] has rigorously shown that the 
groundstate of the ordered Holstein model has to be itin- 
erant, thereby implying that the groundstate can be only 
trapped by an external process, e.g., driven by material 
imperfections. Although it is not clear whether Lowen's 
proof also applies to excited states of the polaron, it is 
generally assumed that all states in the ordered Holstein 
model are itinerant. Hence, if trapping occurs, it has to 
be externally driven. Indeed, translational invariance of 
the Holstein model seems to force all (eigen) states to be 
itinerant. On the other hand, the non-existence of trans- 
lational invariance does not necessarily prohibit itiner- 
ant states. For a single electron in a random potential, 
for example, itinerant states exist, separated from local- 
ized states by a "mobility edge", despite the absence of 
translational invariance. Likewise, as noted above, trans- 
lational invariance per se might be no guarantee that 
all polaron states are itinerant. That is, the important 
question about the trapping properties of excited polaron 
states is essentially unanswered. 

Evidently, the broad theoretical task is to construct 
a rigorous "mobility-edge-theory" for excited polaron 
states, treating extrinsic and (possibly) intrinsic trapping 
processes on an equal footing. In view of the complexity 
of this problem, we are not yet able to develop such a 
rigorous theory. Instead, we present in this paper, as a 
first step, an approximate investigation of the spectral 
and diffusion properties of that part of the spectrum, 
where we expect a rigorous theory to predict a mobility 
edge. In particular, we investigate, in the intermediate 
electron-phonon coupling regime and with an approxi- 
mate technique, the states comprising the high-energy 
edge of the lowest polaronic subband. The dispersion 
is extremely flat in this part of the spectrum indicat- 
ing, perhaps, already the tendency of (intrinsic) trapping. 
Definitely, very small amounts of disorder suffice to drive 
these states into the trapped regime. 

Our investigation is based on the Holstein model 
augmented by site-diagonal, binary-alloy type disorder. 
Measuring energy in units of 2J, where J denotes 
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the electronic (nearest-neighbor) hopping integral, the 
Hamiltonian for the "Holstein alloy" reads in the single- 
electron sector 

i <ij> 

+ S^Xi\i X i\ 

i 

+ nJ2blb,-9Y.(b, + bl)\z><z\. (1) 

i i 

Here, \i > stands for the Wannier state on site i and 
bi[bl] annihilates [creates] an optical phonon on site i. 
The phonon frequency and the electron-phonon coupling 
constant (in units of 2J) are denoted by and g, re- 
spectively. Depending on the occupancy of the site, the 
on-site energy (in units of 2 J) is either or e^. The 
scattering strength S = — can be taken as a mea- 
sure of the disorder. The independent random variables 
Xi are drawn from a bi-modal probability distribution 
p{xi) = cS{xi) + (1 — c)S{xi — 1), with < c < 1 (c 
denotes the concentration of the B sites). 

Spectral properties of the model, such as the density of 
states, the self-energy, and the spectral function, can be 
deduced from the electronic two-point function, which at 
temperature T = and in Wannier representation reads 

g,,{z)xt\G{z)\j >, (2) 

with an averaged one-resolvent 

G(z) = (((0|G(z)|0))). (3) 

Here, G{z) = l/{z — H) denotes the one-resolvent and 
((...)) and (0|...|0) stand, respectively, for averages over 
the quenched disorder and the phonon vacuum. Hence- 
forth, the term "averaged" implies, if not otherwise spec- 
ified, averaged over both quenched disorder and the 
phonon vacuum. 

From the spectral properties alone it is, in general, not 
possible to decide whether the states involved are itin- 
erant (extended) or trapped (locahzed). A more rehable 
method to discriminate between the two kinds of states 
is to investigate the return probability P, i.e., the proba- 
bility to find (in the long time limit) the electron on the 
same site at which it was initially injected. This crite- 
rion has been used, for example, to characterize electronic 
states in a random potential . At T = and in Wan- 
nier representation, the correlation function, whose long 
time limit is the return probability, reads 

P{t)^m\<^W"%><^\e-^"'\^>\m■ (4) 

Using Abel's theorem, P can be conveniently expressed 
in terms of p{2rf), the Laplace transform of l/2P(t/2) 
0. Explicitly, 

P ^ lim P{t) 



= lim 2r]p(2r]) 

r|—^0 

/oc 
— /(uj-iri^uj + iri), (5) 
-oc 27r 

where we have defined a four-point function 

/(zi, za) =< i\K{zi, \i >< i\,Z2)\i >, (6) 

given in terms of an averaged two-resolvent, which (for 
an arbitrary electronic operator Oe) reads 

K{z,,0,,z2) = (((0|G(zi)OeG(z2)|0))). (7) 

The behavior of the spectrally resolved return proba- 
bility, 

P{uj,ri) ^ 2rif{bj ~ iTj.uj + irf), (8) 

in the limit 77 ^ allows one to distinguish trapped 
from itinerant states. In particular, for a trapped po- 
laronic defect state, lim^^o P{u;,r]) is finite, whereas for 
an itinerant polaronic quasi-particle state lim^^o P{i^, v) 
vanishes. 

The exact calculation of the averaged one- and two- 
resolvents for the Holstein alloy is a rather formidable 
task, and, due to lack of rigorous mathematical methods, 
we are forced to adopt an effective single-site averaging 
procedure. Specifically, we employ a generalized dynam- 
ical coherent potential approximation (DCPA) [p"5|-p0[, 
which is known to capture the essentials of the polaron 
formation process. 

Unfortunately, within the DCPA, trapped polaron 
states cannot be unambiguously identified, because 
Imiri^o P{'-^tV) vanishes for aZZ energies uj. Thus, rigor- 
ously speaking, the DCPA predicts all polaron states to 
be itinerant. Nevertheless, it is possible to give a qual- 
itative discussion of the trapping issue. In particular, 
an asymptotic analysis of P{uj,ri) for finite rj allows to 
define a spectrally resolved delocalization time, which is 
the characteristic time scale on which, at a given energy, 
the electron leaves a given site. According to their delo- 
calization times, it is then possibe to distinguish between 
"fast" quasi-particle-like and "sluggish", i.e., temporarily 
trapped, defect-like polaron states. 

The organization of the rest of the paper is as follows. 
Section II contains a general description of the DCPA 
with details about the calculation of the two- and the 
four-point functions defined above. Section III special- 
izes the DCPA formalism to the Bethe lattice with in- 
finite coordination number and discusses representative 
numerical results for the spectral and diffusion proper- 
ties, respectively. Mathematical details not essential for 
the structure of the paper are relegated to two appen- 
dices. Finally, we conclude in section IV with a summary 
of key results and open questions. 
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II. DYNAMICAL COHERENT POTENTIAL 
APPROXIMATION 

A. Perspective 

In contrast to a classical alloy, where an electron ex- 
periences only elastic impurity scattering, the electron in 
a Holstein alloy is subject to elastic impurity and inelas- 
tic phonon scattering. Both scattering processes might 
possibly trap the electron. It is therefore necessary to 
treat them simultaneously. In particular, averages over 
the phonon vacuum and the alloy configurations need to 
be performed on an equal footing. 

A powerful tool to approximately perform these two 
averages is the effective single-site approximation, where 
one embeds an individual scatterer into an effective 
medium to be determined self-consistently by enforcing 
the scattering from a single scatterer to vanish on the 
average. The main advantage of this scheme is the non- 
perturbative treatment of on-site correlations. In partic- 
ular, the non-perturbative treatment of on-site electron- 
phonon correlations is an essential in order to capture po- 
laron formation. The main drawback, on the other hand, 
is the insufficient treatment of inter-site correlations. It 
is the latter which eventually makes lim^_,o 77) to 
vanish for all energies w, irrespective of the model pa- 
rameters, and, thus, prevents a rigorous calculation of 
diffusion and trapping properties. However, as indicated 
in the Introduction, an asymptotic analysis of P{uj, rf) for 
finite 77 partly compensates for this shortcoming. 

In the case of elastic impurity scattering the effective 
single-site approximation yields the well-known coher- 
ent potential approximation (CPA) The DCPA, 
originally developed by Sumi jisj, is a direct extension 
of the CPA to problems involving inelastic scattering. It 
has been primarily used to investigate linear p^ , p^ and 
nonlinear [^-^ optical properties of polaron-excitons 
in molecular crystals. However, it can be applied to any 
on-site inelastic scattering process. Paquet and Leroux- 
Hugon employed the DCPA, e.g., to investigate a 
single quantum particle subject to an on-site potential 
which ffuctuates in time according to discrete or contin- 
uous Markov processes. 

The CPA as well as the DCPA are mean field ap- 
proximations for a single electron subject to impurity 
or phonon scattering, respectively. The most general 
mean field approximation, applicable also to dense quan- 
tum systems with mutual inelastic scattering between its 
constituents, is the dynamical mean field approximation 
(DMFA) H-Hl. [Naturally, in the limit of a single elec- 
tron, interacting with a bath (of phonons or random scat- 
terers), the DMFA reduces to the (D)CPA.] It can be 
shown, assuming that the nearest-neighbor hopping in- 
tegral is scaled by a factor \/\/Z, where Z denotes the 
coordination number of the lattice, that the DMFA pro- 
vides an exact solution of the original model in the limit 
of Z ^ 00. Accordingly, the DMFA [i.e., in the respec- 



tive limits, the (D)CPA] can be either considered as an 
exact theory for lattices with Z — 00 ov as an approxi- 
mate theory for lattices with finite Z. 

The DMFA has been employed by Ciuchi et al. to 
study, in various electron-phonon coupling regimes, the 
groundstate and the spectral properties of a single elec- 
tron in the ordered Holstein model. As indicated above, 
for a single electron, the DMFA reduces to the DCPA. 
The investigations of Sumi and Ciuchi et al. are therefore 
closely related. 

Postponing to section III a discussion of how our in- 
vestigation complements the DCPA work by Sumi 
and the DMFA work by Ciuchi et al. [|ll|, we shall now 
present a derivation of the generalized DCPA formalism 
as needed to determine two-point and four-point func- 
tions for the Holstein alloy. 



B. Calculation of G{z) 

In the spirit of the CPA we introduce an effective 
medium described by an effective one-resolvent 



Geffiz) 
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Heffiz) 



with 



and 



iJe//(z)=i/B+S(z), 



HB = eBj2\i><t\-lYl I* X ■ 



(9) 



(10) 



(11) 



The unknown energy-dependent coherent potential 
(self-energy) operator S(z) defined with respect to the 
reference Hamiltonian Hb will be specified by enforcing 
a self-consistency condition, namely the averaged one- 
resolvent is forced to be equal to the averaged effective 
one-resolvent, i.e.. 



G(z) = (0|G,//(z)|0), 



(12) 



where we have used the configuration independence of 

Geffiz). 

In the above equations, the energy variable z denotes 
the total energy of the coupled electron-phonon system 
and not just the electron energy. The coherent potential 
seen by the electron, on the other hand, depends of course 
only on the electron energy. Therefore, the energy argu- 
ment of the coherent potential is the total energy z minus 
the lattice energy, which is, neglecting the zero point mo- 
tion of the phonons, just the total number of phonons in 
the system times the phonon frequency. A convenient 
way to ensure that the coherent potential is only a func- 
tion of the electron energy is to define, following ref. [po[. 
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a set of projection operators Pq, which project onto the 
phonon subspace with a total number of q phonons, and 
to write for the coherent potential operator 



q i 

= EEV''wi*><*i' 



(13) 



with w*^^) (z) = v{z — qn). From now on we shall adopt the 
convention that any operator or function with a super- 
script (q) has to be taken at the energy z — qft. Note, in 
accordance with the spirit of the CPA, that the coherent 
potential v{z) does not contain any spatial information: 
It neither depends on the site nor on the spatial distri- 
bution of the excited phonons. 

Combining this ansatz for with Eq. (|9|), we obtain 
for the effective one-resolvent 



with auxiliary (electronic) operators 
<?(^)(z)= ^ 



(14) 



(15) 



Apparently, in contrast to the full one-resolvent G{z), 
the effective onc-rcsolvent Geff{z) is a simple operator 
in the phonon Hilbert space. For an arbitrary phonon 
state |{n;}) 



Ge//(z)|M)^ff(S."')(z)|{n,}) 



(16) 



holds, i.e., (phonon) subspaces with different total 
phonon number are decoupled and, as a consequence, the 
self-consistency condition Eq. (^2|) becomes 



G(z)=gW(z). 



(17) 



In order to derive a functional equation for the un- 
known function v{z), which in turn specifies G{z), we per- 
form a multiple-scattering analysis in the product Hilbert 
space of the coupled electron-phonon system. Except for 
the underlying Hilbert space, the mathematical manip- 
ulations parallel the CPA multiple-scattering analysis of 
ref. . The basic relation between the full and effective 
one-resolvents, 

G{Z) = Ge//(Z) + Ge//(z)r(z)Ge//(z), (18) 

involves the total T-matrix, 

T(z) = ^g,(z), (19) 

i 

which, according to ref. |^^, can be expressed in terms 
of individual single-site contributions 



i{z) = ti{z) 



l+Ge//(z)^g,(2 



with the atomic T-matrix given by 

U(z) = [1 - ^H,{z)Geff{z)]-^ AH,{z). 



(20) 



(21) 



From Eq. ( pO| ) we see that Qi{z) is a product of the 
atomic T-matrices ti{z) and an effective wave factor 

^ + Geff{z)J2j^iQji^) describing an effective wave 
incident on site i modified by multiple-scattering events. 
The single-site perturbation 



AH,{z) 



Sxi 



X \i X i\ 



(22) 



comprises bi-modal alloy type fluctuations and coupling 
to localized phonons on site i. 

Combining Eq. ^ with Eq.(|9l) and taking E g. jl^ ) 
into account, the self-consistency condition Eq. (|17|) is 
satisfied if we enforce 



(((0|Q.(z)|0)))=0. 



(23) 



Iteration of Eq. ( p3| ) gives rise to an intractable infinite 
series involving all sites of the lattice. The standard pro- 
cedure to obtain a numerically feasible description is to 
factorize the average in Eq. (|2^) into a product of av- 
erages, i.e., to average atomic T-matrices and effective 
wave factors separately. Relegating a discussion of the 
factorization procedure to the next subsection, we ap- 
proximate Eq. (^ ) by 

(((0|Q.(z)|0))) « (((0|t,(z)|0))) 

X (((0|l + G,//(z)^Q,(z)|0))), (24) 

reducing thereby the self-consistency condition to 
{{{0\ti{z)\0))) = 0. Because phonons on sites j ^ i are 
not affected by scattering events encoded in ti{z) (they 
are "frozen" and act only as "spectators"), the single-site 
self-consistency condition can be formulated more gener- 
ally as 

(((iV,0,|t.(z)|0„iV))) ^ (((0,0.|4^\z)|0„0))) 

= V iV, (25) 

where \ni,N) denotes an arbitrary phonon state with n 
phonons on site i and a total number of N phonons on all 
the other sites. Note the usage of the supersript notation 
defined above. The physical interpretation of Eq. (25) is 
the following: The electron "remembers" how many vir- 
tual phonons it has excited in previous scattering events 
(and, accordingly, how much energy is stored in the lat- 
tice); it does not "remember", however, at which sites 
these events took place. 
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Instead of working directly with Eq. (25), it is more 
convenient to introduce an equivalent polaron-impurity 
model (PIM) which describes a single perturba- 

tion AHi{z) self-consistently embedded into an effective 
medium. Mathematically, the PIM is formulated by 

D,{z) = G,ff{z) + Geff{z)U{z)Geff{z) (26) 
= G.ffiz) + G,ffiz)AH,{z)D,iz). (27) 



The self-consistency equation ( pSD is then transformed 
into 

(((0,0.|i?f\z)|0„0)»=ffW(z), (28) 
which in turn yields for the averaged one-resolvent 

G{z) = g(°\z) = {{{OMmmM)- (29) 

As shown in appendix A, it is straighforward to derive 
from the PIM a non-linear functional equation for the 
local two-point function. We get 



Qiiiz) 



1 - c 



[F^'(°)(2)]-i 



[F-^.W(z)]-i- 



2g^ 



(30) 



[F^^(i)(z)]-i 



25' 



with 



z — — e\ + Gii^iz — Tifi) — TZ[Qii{z — nfl)] ' 



(31) 



where TZ[$] denotes the inverse Hilbert transform (see 
below) . 

For c = Eqs. (|^)-(|l|) have been given before by 
Sumi ||l^ and independently by Ciuchi et al. [|ll|. Various 
limits of the c = equations, e.g., weak, intermediate, 
and strong electron-phonon coupling limits, have been 
discussed in these references. The non-linear functional 
couples local two-point functions with energies shifted by 
any negative integer multiple of f2. This accounts for the 
fact that at zero temperature only phonon emission is 
possible. At finite temperature, in contrast, absorption 
of thermally excited phonons is also allowed. In that 
case, local two-point functions with energies shifted by 
any integer multiple of Q, negative as well as positive, 
would be coupled | pT| , p^ . 

The coherent potential v{z) can be directly obtained 
from Eq. (|29|). Specifically, introducing the Hilbert 
transform Ti.[^] corresponding to the bare density of 
states iVo(£;) = {1/N) J2kS{E - e^:), the local two-point 
function can be written as 



No{E) 



z - EB ~ v(o) (z) - E 
= n[z-eB-v^"Hz)]. (32) 

Employing the inverse Hilbert transform, defined by 7?. * 
Ti-l^] = finally leads to 



(33) 



The local two-point function Gu^z) uniquely determines 
the coherent potential v^^\z). Once w'^''^(z) is known 
any matrix element of the averaged one-resolvent can be 
obtained. 



C. Validity of the Factorization Procedure 

In the previous subsection we employed a factorization 
procedure to reduce the multi-site self-consistency condi- 
tion to a single-site condition. Clearly, the mathemat- 
ically exaxt average of the operator Qi contains many 
terms which are not zero, even if the atomic T-matrix 
vanishes. The single-site self-consistency condition can 
be only approximate. We now analyze the validity of the 
factorization procedure taking advantage of the connec- 
tion between the DCPA and the DMFA. As a by-product, 
we shall find an efficient way to derive the DCPA equa- 
tions for the averaged two-resolvent. 

To that end, we follow ref. and rearrange the 
multiple-scattering series into clusters containing a fixed 
number p of lattice sites. We write (suppressing the en- 
ergy variable z) 



(34) 



with Q^^ — ti and 

H Gefft^^_^{iJl,...,jp-2) 



(35) 



forp > 2. The cluster T-matrices t^^'_^(z, ji, jp_2) are 
defined as the sum of all scattering processes involving 
at most p different sites i,ji, jp-i with the constraint 
that the entrance site has to be jp-i- This constraint 
ensures that all p sites contributing to q'^' are connected. 
The exit site, on the other hand, can be any site of the p- 
cluster. Notice that, in contrast to the original multiple- 
scattering series given in Eq. (|2^), where only successive 
summation indices had to be different, all summation 



indices in Eq. (35) are different. 

We are now in a position to evaluate the importance 
of each individual multiple-scattering term contributing 
to Qi. Keeping in mind that all processes where a site 
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j is visited only once vanish due to the single-site self- 
consistency condition Eq. (p5|), only processes where all 
sites are at least visited twice require further analysis. 
Figure [I] schematically depicts the first non- vanishing p = 
2 process (i.e., a fourth order process in which each of the 
two sites is visited twice). Writing 



=<z|(0,n,|irV.,0)|z> 



Wi 



(36) 



for the (configuration dependent) matrix elements of the 
atomic T-matrix, this process is analytically given by 



Q 



i — ''Or ^ij '•Oq ^ji ''rO ^ ij ''gO 



(37) 



To proceed further, we now imagine at this point that 
the hopping matrix element is scaled by l/^/Z (as in the 
case of the DMFA). Then, we immediately see that for 
large Z Eq. ( |37| ) is at least ~ 1/ \fz because it contains 
three off-diagonal two-point functions and only one free 
summation over sites j, i.e., for Z ^ oo this process does 
not contribute to Qi at all. In fact, it can be shown that 
all p-cluster processes in which all p sites are at least 
visited twice are suppressed by a factor containing the 
inverse of the coordination number Z. As a consequence, 
for Z — > cxD these processes vanish and need not be ex- 
plicitly considered in the averaging procedure. 

The implication of the above discussion is three-fold. 
First, it corroborates, using multiple-scattering theory, 
the conventional CPA technique, previous work by Vlam- 
ing and VoUhardt |^^, who showed, acknowledging ear- 
lier work by Schwartz and Siggia |2^, that for a lattice 
with infinite coordination number any CPA-type theory 
becomes exact, i.e., for infinite coordination number, the 
single-site self-consistency condition is indeed sufficient 
to ensure that the average of Qi as a whole vanishes. 

Second, it indicates that even for a finite coordination 
number, processes which do not vanish due to the single- 
site self-consistency condition are at least suppressed. 

Third, from a formal point of view, the DCPA is equiv- 
alent to the replacement of the full operator Qi by a re- 
duced operator Qf^^ containing only self-avoiding paths 
(SAP); Qf^^ can be obtained from Eq. ^ by set- 
ting t^j__^{i,ji, ...,jp^2) — tjp-^ P- The coherent poten- 
tial v{z) is then adjusted in such a way that the aver- 
aged Qi^^ vanishes. This procedure can be envisaged, 
in physical terms, as a loss of spatial memory^ because 
the electron "assumes" that each site it encounters is vis- 
ited the first time, i.e., the electron does not "remember" 
whether it has visited the site before. Naturally, for a 
lattice with infinite coordination number, this is exact, 
because the probability that the electron returns to the 
same lattice site is vanishingly small. 

The replacement Qi Qf^^, i.e., the erasur of the 
spatial memory of the electron, is of course consistent 
with the ansatz for the (local) self-energy operator. On 
the other hand, this replacement leads to an insufficient 
treatment of inter-site processes in the calculation of the 



averaged two-resolvent and, eventually, prevents a rigor- 
ous calculation of the diffusion and trapping properties. 

Keeping this shortcoming in mind, we shall neverthe- 
less employ this formal replacement from the outset in 
the next subsection to simplify the calculation of the four- 
point function. 



D. Calculation of K{zi,Oe, Z2) 

The discussion of the previous subsection suggests that 
the single-site self-consistency condition is, from a formal 
point of view, equivalent to the replacement Qi — > Qf^^, 
i.e., only SAPs need to be explicitly considered in the 
averaging procedure. We now show that this replacement 
allows a rather straightforward derivation of the averaged 
two- resolvent. 

For that purpose, we first recast Eq. (^, using Eqs. 
( |l6| ) and (fsj) together with the self-consistency condi- 
tion, Eq. (pq), into 

K(zi,0e,Z2) = .g("H^l) [Oe+fizi,Oe,Z2)] 

X .g(")(z2). (38) 

Keeping only site-diagonal terms, the electronic vertex 
operator is given by 

f (Zi, Oe, Z2) = ^(((0|r,(zi, Oe, ^2)|0))), (39) 
i 

with a local vertex operator 



F i{zi, Oe, Z2) = 

Q,^'4^(zi)Ge//(zi)OeGe//(z2)QP''(^2). 



(40) 



In the above equations we have kept only retrace- 
able SAP's where the second leg of the path, encoded 
in Qf^^, retraces in reversed order the first leg of the 
path, encoded in Qf^^ and neglected retraceable SAP's 
where the second leg retraces the first leg in the same 
order. In more familiar terms (see Fig. we have kept 
the ladder diagrams and neglected (maximally) crossed 
diagrams. This choice is consistent with the coherent 
potential determined in the previous subsection. In par- 
ticular, it can be shown that the local vertex defined in 
the ladder approximation and the local self-energy obey 
the Ward identity corresponding to particle conservation 

il- 

In analogy with the classical alloy problem, we expect 
the maximally crossed diagrams to yield a significant (for 
certain model parameters possibly diverging) contribu- 
tion to the local vertex function. As a consequence, a 
consistent treatment of these processes might render the 
potential to force lim^_,o P(w, rj) to be finite for some en- 
ergies uj and for certain model parameters. Accordingly, 
itinerant and trapped polaron states could be rigorously 
distinguished. We did not, however, succeed in taking 
them into account. 
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To derive an explicit equation for the focal vertex func- 
tion (in the ladder approximation) , we have to perform 
the average of Eq. (^0|). To that end, it is advantageous 
to derive a formal operator equation for the local ver- 
tex operator F^. For simplicity, we shall from now on 
suppress the arguments of the various operators. The 
implied arguments should be clear from the context. If 
we visualize Eq. (|40| ) in terms of diagrams (see Fig. ||), 
we readily obtain 



ff 



GeffU, 



(41) 



vertex functions with energies shifted by any negative in- 
teger multiple of Vl. 

Employing spatial Fourier transforms. 



Q ^ 



(r) 



('■) 



-iQRi 



< 



4' - E e-^'5-[«.-K.l^W^(0 



E' 



(46) 
(47) 
(48) 

(49) 



i — i) 

where F^- stands for the vertex operator on site j given 
by summing up all "ladder-type paths" not visiting site 
i. Consequently, as far as the configuration average is 
concerned, the atomic T-matrices ti and the term be- 
tween the brackets are statistically independent and the 
configuration average factorizes, i.e.. 



((r,:)) = ((i,Ge//(( 



))Ge//t.)). (42) 



For a macroscopic solid we expect ((F^- essentially 
to be the same operator as {(Tj))- Hence, if we dealt 
only with alloy-type disorder, Eq. (|2|) would be the 
result given by Velicky |^^. To cope with inelastic scat- 
tering, we have to keep in mind, however, that the vertex 
function (similar to the coherent potential) depends only 
on the electron energy, i.e., we have to keep track of the 
energy stored in the lattice (given by the total number 
of virtual phonons). Consistent with the ansatz for the 
local self-energy operator, we therefore write 



((r.»-EAr, 



(9) 



(43) 



with j^'^^ = lii^i ~ <l^,Oe,Z2 — q^)- Thus, the phonon 
vacuum expectation value of Eq. (E3) becomes 



7^ 



(0) 



((E^^°^ 



Oq 



K. 



(9) 



Q{q) {i)Q{q) 



with 



(44) 



(45) 



Note, the only configuration dependent quantities are 
now the matrix elements of the atomic T-matrix. 

The linear functional for the vertex function is not "lo- 
cal" in the energy variables zi and Z2- Instead, due to 
the possibility of phonon emission, Eq. (44) couples all 



it is straigthforward to combine Eqs. (^4|) and (^5|) and 
to derive, after some rearrangements, a single equation 
for the Fourier-transformed local vertex function: 



(9) 

Q 



(9) 

Q 



(50) 



This equation reduces in the limit = to Velicky's 
vertex equation ||2^ for a model with a Gaussian distri- 
bution of on-site energies on top of a bi-modal distribu- 
tion of on-site energies. Note, moreover, that 7^*-' may 

be used to calculate any matrix element of the averaged 
two-resolvent of interest, i.e., the averaged two-resolvent 



(0) 



is completely determined by 7 



III. NUMERICAL RESULTS FOR THE Z = oq 
BETHE LATTICE 

A. Final Set of Equations 

In the previous section, we neither specified the co- 
ordination number Z nor the structure of the underly- 
ing lattice. The DCPA equations for the one- and two- 
resolvents can be used for arbitrary lattices. On the other 
hand, we have shown in subsection II. C that the single- 
site factorization procedure essentially ignores processes, 
where the electron returns to the same lattice site. It 
is therefore consistent with the DCPA to use a lattice, 
which, by construction, completely prohibits this type of 
processes. A particularly simple case is the Z = 00 Bethc 
lattice. (It should be emphasized, however, that, within 
the DCPA, the numerical results, in particular for the 
spectrally resolved return probability, are not affected 
much by the chosen lattice.) 

The main simplification, as far as the calculation of 
the local two-point function is concerned, arises from the 
fact that the inverse Hilbert transform for a Z = 00 Bethe 
lattice is an algebraic relation, namely TV\^ ~ 1/4^-1-1/^, 
which makes the self-consistent solution of Eqs. (30)-(|3l|) 
particularly easy. 
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The calculation of the local vertex function simplifies 
even more, because, for a Z = oo lattice, Eq. (Hi) re- 
duces (for almost all Q) to ^^^(2:1, Z2) ~ Qu\zi)g^ {Z2) 

p6| . Therefore, the second term on the rhs of Eq. (|5C 
vanishes and the local vertex function becomes 



(0) 



(«) 



(51) 



which yields for the Fourier transform of Eq. ([45[) 



K^=o^+G^'G^'j:i(fX))<^- (52) 



To investigate the trapping properties of the polaron 
states, we need the four-point function /(zi,Z2) defined 
in Eq. (||). To that end, we set Og — |* >< i\ in the 
above equation and obtain after Fourier transformation 

f{zi, Z2 ) = Qvi{zi)Gvi{z2) + Qii{zi)Qvi{z2) 

X Y.^(f^^^Z^^{z2)))g'^{z^g'^{z2). (53) 

q 

This equation holds for any lattice with Z — 00. The 
structure of the lattice, i.e., in our case, the structure of 
the Bethe lattice, enters only through the local two-point 
function Qii{z). 

In appendix B we show that Eq. (^3|) can be brought 
into a form which allows us to calculate /(w — irj^ui -\- irj) 

from the amplitudes V^q^^ (uj + irj) defined in appendix 
A. Specifically, 

oc 

f{u; -^rJ,u; + ii^) = (1 - c) ^ \vff^ {u; + 

00 

+ cJ2\T^T\^ + ^v)\', (54) 

i.e., the calculation of the four-point function /{uj—irj, UJ+ 
irj) requires the same (moderate) numerical effort as the 
calculation of the local two-point function Qu {u + irj) . 

The numerical strategy is now as follows. First, we 
truncate the continued fraction expansions in Eq. ( |30| ) 

and find the fixed point Q\1^^\uj + irf) by iteration. The 
density of states then follows directly from 



N{uj) 



(55) 



To study the individual contributions of the A and B 
sites, we also calculate the component density of states 



(56) 



The amplitudes 'D'^^'^^ , which are also needed in the cal- 
culation of the spectrally resolved return probability (see 
below), are given by the first (A = A) and the second 



(A = B) term on the rhs of Eq. (|3q ) without the factors 
1 — c and c, respectively, and with Qii{z) = Q\1^^\z). 

Second, we insert q\(^^\ijJ + in Eq. (|33| ) and obtain 
the coherent potential v^^^u -f zry), which yields, using 
Eqs. ([isl ) and (17), the spectral function 



- —Ini 

TT LU + irj — eg 



1 



v^^") {uj + irj) 



(57) 



Third, to calculate the spectrally resolved return prob- 
ability PU£, 1]) defined in Eq. (^) , we truncate the sums 
in Eq. ( p4[ ) and successively construct the amplitudes 
V^Q^^uj + irj) from the recursion relation, Eq. ( [A1C| ), 

starting from I?o6^°^ (w -|- iry) with X = A, B. The trun- 
cation procedures and the subsequent numerical solution 
of the truncated equations for small rj > converge very 
well. 



B. Spectral properties 

In this subsection we discuss the spectral properties 
of a subclass of polaron states, which we expect to be 
very susceptible to an (intrinsic or extrinsic) trapping 
transition. However, as mentioned in the Introduction, 
spectral properties alone cannot be used, even within a 
rigorous theory, to decide whether states are itinerant or 
trapped. 

We start with the ordered Holstein model [model (|l|) 
with 5 = 0], which is characterized by two parameters, 
A = 2g'^/Q and a = g/V,. (Recall, we measure energy 
in units of 2J.) These parameters are conventionally 
used to define strong, intermediate, and weak electron- 
phonon coupling regimes The overall DCPA spec- 

tral features of which have been thor oug hly investigated 
by Ciuchi et al. ||ll| and by Sumi [p^ . The most in- 
triguing of these regimes is, perhaps, the intermediate 
electron-phonon coupling regime, approximately defined 
by A > 1 and a > 1, where a finite number of polaronic 
subbands starts to develop on the low energy side of the 
bare density of states. The formation of these polaronic 
subbands, in particular as a function of A and a, has 
been intensively studied in above mentioned references. 
For our purpose, it is sufficient to focus, for a fixed A and 
a, on the lowest polaronic subband. 

Generic numerical results for the ordered Holstein 
model in the intermediate electron-phonon coupling 
regime are summarized in Fig. ^. The polaron param- 
eters are given in Table |. We first focus on panels (a) 
and (b), which depict the density of states N{u!), the 
imaginary part of the coherent potential Imv{uj), and 
the spectral function A{uj^e-j^, respectively. The spec- 
trally resolved return probability P(a;, rf) shown in panel 
(c) will be discussed in the next subsection. 

The main spectral characteristics in the intermedi- 
ate electron-phonon coupling regime is the formation 
of a well-established lowest polaronic subband with an 
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asymmetric, steeple- like density of states, separated from 
the rest of the spectrum by a hard gap. Below the 
phonon emission threshold, which is inside the second po- 
laronic subband, no residual scattering takes place, i.e., 
Imv{uj) = 0. That is, within the lowest polaronic sub- 
band, the electron and the phonons arrange themselves 
into a new (composite) entity - a polaronic quasi-particle 
0. [The pole of Imv{uj) slightly above the high-energy 
edge of the lowest polaronic subband does not imply a 
diverging scattering rate; it merely signals the presence 
of the gap.] That the quasi-particle concept indeed ap- 
plies can be also seen from the spectral function, which 
features a (dispersive) quasi-particle peak in the energy 
range of the lowest polaronic subband. Notice the pro- 
nounced band flattening and the vanishing quasi-particle 
weight in the flat part of the dispersion. The band flat- 
tening is the reason for the steeple- like density of states. 
At higher energies a second quasi-particle peak corre- 
sponding to the second polaronic subband is visible until 
the point where it merges with the one-phonon scattering 
continuum. 

Both the steeple-like shape of the density of states and 
the band flattening at the high-energy edge of the lowest 
polaronic subband are an indirect consequence of the pole 
in the imaginary part of the coherent potential. As can 
be seen in Figs. |^ and the pole in the imaginary part 
is accompanied by a divergence in the real part. [Both are 
broadened due to the finite value of r?; the pole becomes 
a Lorentzian and the divergence becomes a pronounced 
kink.] The divergence, in turn, induces a pronounced en- 
ergy dependence of Rev{uj) just below the high-energy 
edge and it is this strong energy dependence which, as 
we shall show below, eventually yields the steeple-like 
density of states and the band flattening. Physically, the 
energy dependence of the real part of the coherent poten- 
tial reflects the increasing phonon admixture in the states 
as the high-energy edge of the subband is approached. 

The discussion of the steeple-like subband density of 
states is facilitated by considering 



Nist{Lo) = -Re^l -[u-CB- Rev{uj)Y, (58) 

TT 

an exact analytical expression for the lowest subband 
density of states, which follows from Eqs. (55) and (^) 
in the limit ry ^ and the fact that for the lowest sub- 
band Imv{uj) = 0. Two conclusions can be drawn from 
Eq. (|5^). First, the position cj,„ of the maximum of 
the subband density of states does not coincide with the 
high-energy edge toh, because they are respectively given 
by uj„i — es — Rev{uJm) — (thin solid line in Fig. @a) 
and uJh — £b — Rcv{ujh) = 1 (thin dashed line in Fig. |5p,). 
Hence, the density of states at the high-energy egde is, 
albeit very step-like, not discontinuous, even for — > 0. 
Second, Eq. ( p8| ) together with Figs, ^a and ^3 reveal 
that the rapid change of Rev(uj) between ujm and w/,,, 
caused by the divergence above ujh, produces the rapidly 
dropping density of states in this particular energy range, 
and hence, the steeple-like shape. 



Figures ||a and ^ suggest, moreover, that the rapidly 
changing Rev{Lo) is also responsible for the band flatten- 
ing. To be more specific, we recall that the dispersion 
^{^j:) is given by the poles of the spectral function de- 
fined in Eq. (M) . Therefore, the change of the dispersion 



with becomes 



(1 



dRev{Lo) 
cLu 



lt^='^(ej)) 



(59) 



which is, of course, particularly small for a;,„ < < 
LOfi, since the divergence yields a large derivative of 
Reviuj). It is, therefore, the divergence-induced strong 
w-dependence of the real part of the coherent potential 
which not only causes the steeple-like subband density of 
states but also the pronounced flattening of the subband 
dispersion. 

The DCPA results for the lowest polaronic subband 
corroborate recently obtained direct numerical simula- 
tion results for a finite Holstein model in the single- 
electron sector 10-^ . Specifically, the band flattening and 
the vanishing quasi-particle weight in the flat part of the 
dispersion are in good qualitative agreement. Both ef- 
fects, characteristic of the intermediate electron-phonon 
coupling regime, are manifestations of the strong on-site 
electron-phonon correlations, which yield an increasing 
phonon admixture in the states comprising the high- 
energy edge of the lowest polaronic subband § . 

Clearly, the phonon admixture in the states must sig- 
nificantly affect the diffusion properties. In particular, 
because the dispersion becomes extremely flat in the 
vicinity of the high-energy edge of the lowest polaronic 
subband, these states might be temporarily trapped. In 
the next subsection, we will present an analysis of the 
spectrally resolved return probability P{uj,ri), which in- 
deed shows that the high-energy edge states are ex- 
tremely "sluggish" . 

However, first we illustrate the modifications of the 
lowest polaronic subband due to alloying. We expect 
disorder on the scale of the subband width to show the 
most dramatic effects and therefore restrict S [in Eq. (|^)] 
to small values. Further, to be specific, we chose the B 
atoms to be energetically close to the high-energy edge 
of the subbands, i.e., we consider 6 < 0. 

Figure ^ shows representative numerical results for the 
Holstein alloy. The polaron parameters are the same as 
in Fig. ^ (see Table |) and the alloy parameters are c = 
0.5 and 6 = —0.01. First, we focus again on panels (a) 
and (b), which summarize the spectral properties in the 
vicinity of the lowest polaronic subband. The spectrally 
resolved return probability shown in panel (c) will be 
discussed in the next subsection. 

The main effect of alloying is the formation of two 
mini-subbands. For the chosen parameters both mini- 
subband density of states exhibit a pronounced asym- 
metric, steeple-like shape (see Fig. ||a). From the strong 
asymmetry of the component density of states, which are 
also depicted in Fig. we find, moreover, that, al- 
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though both A- and B-sites contribute almost equally to 
both mini-subbands, the steeple-like shape of the mini- 
subband density of states comes entirely either from the 
A-sites (lower mini-subband) or from the B-sites (up- 
per mini-subband) . As in the case without disorder, the 
steeple-like shape is due to strong on-site electron-phonon 
correlations, which yield a large phonon admixture in the 
states at the high-energy edges of the respective mini- 
subbands, and, as a consequence, to rather flat disper- 
sions in these parts of the spectrum. The band flattening 
can be clearly seen from the spectral function shown in 
Fig. ^. 

The detailed appearance of the mini-subbands depends 
on the concentration c and the scattering strength 6. To 
indicate the concentration dependence, for example, we 
show in Fig. |^ the density of states N{uj) for a fixed 
scattering strength S = —0.004 and concentrations vary- 
ing from the pure A crystal (c = 0) to the pure B crystal 
limit (c = 1). Note, first, the noticeable distortions of 
the density of states (at the high-energy edge) despite 
the fact that S is much smaller than the polaronic sub- 
band width and, second, the lack of c 1 — c symmetry. 
That rather small amounts of disorder suffice to restruc- 
ture the high-energy edge can be also seen in Fig. ||, 
where the lower panel shows the density of states and 
the imaginary part of the coherent potential for c = 0.5 
and three different scattering strengths 6. 

It is not necessary to present an exhaustive investi- 
gation of the alloying-induced spectral changes in the 
whole c — S parameter plane. The data presented suf- 
fice already to show that the modifications of the lowest 
polaronic subband are very similar to the ones found in 
an AB alloy with an artificial steeple-like bare density of 
states ||3l| . Indeed, the broken c <-> 1 — c symmetry, the 
appearance of a B mini-subband at scattering strengths 
much smaller than the polaronic subband width, and the 
asymmetry between the component density of states have 
their counterparts in this artificial steeple alloy. 

In the artificial steeple alloy, as well as in the Holstein 
alloy, states comprising the steeple-like structure in the 
density of states are very susceptible to impurity scat- 
tering. For the Holstein alloy this is intuitively clear, 
because the steeple in the density of states contains the 
states corresponding to the flat part of the dispersion. 
As indicated above, due to the large phonon admixture, 
these states might be temporarily trapped, which sug- 
gests that they feel the on-site, alloy-type disorder much 
more strongly than the other states. The discussion of 
the diffusion properties in the next subsection will sub- 
stantiate this intuitive picture. 

C. Diffusion properties 

The spectral properties in the vicinity of the high- 
energy edge of the lowest polaronic subband suggest that 
this part of the spectrum is very susceptible to an (in- 



trinsic or extrinsic) trapping transition. 

In the Introduction, we emphasized that the problem 
of trapping of a polaron should be methodologically ad- 
dressed with the same techniques used to investigate the 
problem of localization of an electron in a random poten- 
tial. In particular, we stressed that lim^_»o -P(w, rj) would 
be a rigorous criterion to distinguish itinerant polaronic 
quasi-particle states from (any) trapped polaronic defect 
states. We also pointed out, that a rigorous investigation 
of the trapping issue is unfortunately beyond the DCPA 
because, within the DCPA, lim^^o Pi'^, v) vanishes for 
all energies lo. Although trapped states cannot be unam- 
biguously identified, it is possible to identify temporarily 
trapped states through the 7y-asymptotics of P{Lu,ri). 

To motivate the investigation of the jy-asymptotics of 
P(cj,77), we first emphasize that the spectrally resolved 
return probability P(w, rj) contains valuable information 
about the diffusion properties not only in the limit rj = 
(which is usually considered ||l^) but also for a finite 
77 > 0. For a correct interpretation of the data, we have 
to keep in mind, however, that a finite rj broadens all 
spectral features. In particular, the pole [divergence] in 
Imv{uj) [Rev{uj)] becomes, as mentioned in the previous 
subsection, a Lorentzian [pronounced kink] , smearing out 
the high-energy edges of the respective bands. A finite 
77 yields therefore small artificial tails extending into the 
gap regions. These tails have no direct physical meaning. 
In the following we analyze, therefore, the 77— asymptotics 
of P{w, rj) only for energies lo below the respective sharp 
high-energy edges obtained for 77 ^ (see Fig. ||d and 
the insets of Figs. ||-|n]). 

The physical content of P(lo, 77) can be extracted from 
Figs. |d, |c, and the upper panel of Fig. ||, where 
we plot, for the model parameters given in the respec- 
tive captions, P{uj,ri) for rj — 10^^. All the data show 
that P{uj, r/) does not uniformly (in to) approach zero as 
77 — > 0. Instead, P{uj, rj) is substantially enhanced for 
certain energies oj. Close inspection of the data shows, 
moreover, that the maximum of P{io,r]) always occurs 
slightly below the high-energy edges of the respective 
(mini-)subbands, i.e., in the spectral regions where the 
dispersions become extremely flat. (This can be most 
clearly seen in Fig. |[) The enhancement of P{uj,rf) is 
caused by the band flattening and is therefore a direct 
measure of the sluggishness of the corresponding states 
due to their large phonon admixture. More precisely, 
because I/77 can be interpreted as a characteristic time 
scale, a large value of P(w, 77) implies that the electron 
in the state at energy u) is not yet delocalized on the 
time scale 1/77. Thus, the 7y— asymptotics of P{uj,r]), i.e., 
the way P{uj,ri) approaches zero for 7; — > 0, reveals the 
time scale on which the state at energy lu is temporarily 
trapped. [Note, once more, we study the 77— asymptotics 
of P{u!,r]) only for co below the 77 — + high-energy edge 
and not for energies belonging to the artificial tail which 
is present for a finite 7;.] 

In what follows we put these qualitative considerations 
on a mathematical basis. To that end, we first note that 



10 



Eq. (H) suggests to interpret f {uj — ir] , lj + ir]) as a "spec- 
tral density" corresponding to p{2ri), the Laplace image 
of l/2P(t/2). That is, f{u! — irjjUj + irj) gives the spec- 
trally resolved diffusion behavior of the electron in the r] 
domain. Accordingly, the counterpart of /(w — irj, uj + irf) 
in the time domain, which we denote by F{uj,t), is the 
"spectral density" corresponding to l/2P{t/2). That is, 
it gives the spectrally resolved diffusion behavior of the 
electron in the time domain. 

That the ?7-asymptotics of P{lo, rj) indeed contain the 
spectrally resolved diffusion behavior of the electron can 
be clearly seen in Fig. ||, where, for the ordered Hol- 
stein model with the polaron parameters given in Table 
H P{iLJ, rf) is plotted as a function of rj for three represen- 
tative energies uj within the lowest polaronic subband, 
i.e., for energies below the 77 — > high-energy edge. For 
the Holstein alloy, the overall behavior of P{LO,rf) is the 
same. All three curves feature a hump whose shape and 
position strongly depend on u. For w in the immediate 
vicinity to the high-energy edge of the polaronic subband, 
the hump appears at a very small r/ and is very sharp. 
Far inside the subband, in contrast, the maximum ap- 
pears at a larger rj and is less pronounced. The distinct 
Lo dependence of the hump indicates different charateris- 
tic time scales for the diffusion behavior of the respective 
states. 

To explicitly extract this time scale, we perform an 
asymptotic expansion of 77) for small rj. The data 
shown in Fig. | suggest making the following ansatz 



(1 



(60) 



and determining the parameters ??o(w), and v{lo 

through a least-square fit. The dashed lines in Fig, 
show that this approach works reasonably well. Using 
the definition of 77) given in Eq. (||) and employing a 
general mathematical theorem about Laplace transforms 
p2| , we obtain the asymptotic behavior of F{uj,t) for 
large t [note, since we measure energy in units of 2J, 
time is mearured in units of (2J)^^]: 



f,„.0»A-(.)||2M)^.-.-"., 

2r(l + v(yj)) 



(61) 



where T{x) denotes the Gamma function. 

From this equation we infer that the parameter »7o(t^) 
gives rise to a characteristic time scale. Defining T^{lo) = 
1/770 (w), we see, in particular, that for t ^ Tq(ljj) 
the function F{uj,t) decays exponentially to zero. For 
t <^ Tq{uj), however, F{uj,t) is finite, implying that the 
state at energy uj is not yet delocalized. It is therefore 
natural to interprete Tq{ijj) as a (spectrally resolved) de- 
localization time characterizing temporary trapping, i.e., 
Tq{u!) quantifies the "sluggishness" of the state at energy 

UJ. 

In Fig. ^we plot To{uj) for the ordered Holstein model 
in the vicinity of the high-energy edge of the lowest po- 
laronic subband. The polaron parameters are the same 



as in Fig. ^ (see Table |). In the immediate vicinity 
of the high-energy edge, where the dispersion becomes 
extremely flat due to the high phonon admixture in the 
states, the delocalization times increase roughly three or- 
ders of magnitude (but do not diverge). That is, the 
states in the flat part of the polaronic subband appear to 
be temporarily trapped on the time scales where the rest 
of the states are already delocalized. Clearly, the band 
flattening and the enhanced delocalization times are cor- 
related and have, moreover, the same microscopic origin, 
namely the increased phonon admixture in the states at 
the high-energy edge of the subband. More formally, it is 
the divergence-induced strong oj— dependence of Rev{uj) 
for ujm < < (^/i which causes the large delocalization 
times at the high-energy edge of the lowest polaronic sub- 
band. 

At this point it is appropriate to comment on the work 
of Hotta and Takada [|3| who suggested that for energies 
where the self-energy (or, in our case, the coherent po- 
tential) diverges, states would be trapped because of the 
infinitely strong electron-phonon coupling. They stud- 
ied the half-filled (ordered) Holstein-Hubbard model, but 
their conclusion is quite general and should, if valid, hold 
for any model which produces divergences in the self- 
energy. Such divergences are, however, always inside a 
gap. In fact, they signal the presence of a gap, as can be 
seen, e.g., in Figs. || and |[ that is, there are no states at 
the divergence energy which could be trapped. In con- 
trast, our investigation suggests that, it is the effect the 
divergence has on states within the spectrum, namely, 
the band flattening it causes below the high-energy edge 
of the lowest polaronic subband, which could perhaps in- 
duce trapping in that particular energy range. Within 
the DCPA, however, the effect is too weak, giving rise 
only to large but finite delocalization times and, hence, 
only to temporary trapping. 

It has been also suggested that, if, due to some 
mechanism, a finite density of states appears at the diver- 
gence energy, the corresponding states would be trapped. 
A finite rj simulates such a mechanism. However, as can 
be seen in Fig. ^ the maximum of P{uj, rj) does not occur 
in the tail but roughly at the maximum of the subband 
density of states. Therefore, at least within the DCPA, it 
is the band flattening and not the tail which is responsible 
for large delocalization times. 

The concept of the delocalization time provides an ap- 
pealing explanation of the spectral changes induced by 
alloying. In the previous subsection we have seen that 
the high-energy edge of the polaronic subband responds 
much more strongly to alloying than the bottom of the 
subband. From the perspective of the delocalization time 
this is now clear, because the high-energy edge contains 
"sluggish" defect-like states, whose large delocalization 
times make them feel the randomness of the on-site po- 
tential particularly strongly. 

The delocalization times are also modified due to al- 
loying. This is shown in Fig. ^ for a particular set of 
alloy parameters, c = 0.5 and S — —0.002. The polaron 
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parameters are the same as before (see Table |). In the 
previous subsection we have seen that the main effect 
of alloying is the formation of two mini-subbands with, 
due to the large phonon admixture in the states, rather 
flat dispersions at the two high-energy edges. As a con- 
sequence, the delocalization times are now enhanced in 
two narrow energy regions: At the high-energy edge of 
the A mini-subband and at the high-energy edge of the 
B mini-subband. Note, as in the case without disorder, 
the enhancement of the delocalization times is due to the 
band flattening effect, i.e., it is driven by strong on-site 
electron-phonon correlations. Even with the assistance 
of electron-impurity scattering, the on-site correlations 
are however not strong enough to produce divergent de- 
localization times. 

The numerical results for the delocalization times sug- 
gest that we may tentatively distinguish two types of 
states: "Fast" quasi-partice-like states and "sluggish" 
defect-like states. In the ordered Holstein model, defect- 
like states appear at the high-energy edge of the pola- 
ronic subband. Whereas in the Holstein alloy, defect-like 
states occur at the high-energy edges of the two mini- 
subbands. With and without disorder, strong on-site 
electron-phonon correlations are responsible for the for- 
mation of defect-like states. Because the delocalization 
times of the defect-like states are several orders of mag- 
nitude larger than the delocalization times of the quasi- 
particle-like states, it is moreover conceivable that the 
sluggishness of the defect-like states may anticipate the 
onset of trapping. 

IV. CONCLUSIONS 

We have provided a self-contained description of the 
DCPA and employed it to investigate the dynamics of 
a single electron in the Holstein model augmented by 
site-diagonal, binary-alloy type disorder (Holstein alloy) . 
Using multiple-scattering theory, the conventional CPA 
technique, we derived the DCPA equations for the aver- 
aged two- and four-point functions and emphasized the 
main approximation involved: namely, after each hop 
from one lattice site to another the electron's spatial 
memory is erased. We visualized the lack of spatial mem- 
ory in terms of self-avoiding paths and employed this 
picture to give an intuitive explanation why the DCPA 
becomes exact for a system with infinite coordination 
number. 

Our numerical results focused on the intermediate 
electron-phonon coupling regime, where polaronic sub- 
bands start to emerge. We investigated, for representa- 
tive parameter sets, the effect of alloying on the high- 
energy edge of the lowest polaronic subband. These 
states are very susceptible to impurity scattering because 
of their large phonon admixture, which leads to a flat 
dispersion in this part of the spectrum, and, as a conse- 
quence, to a stecple-like density of states. The observed 



modifications of the spectral properties are therefore rem- 
iniscent of typical CPA results obtained for a model with 
a steeple-like density of states: A lack of c ^ 1— c symme- 
try, an appearance of mini-subbands in the vicinity of the 
steeple (i.e., high-energy edge) of the subband for scat- 
tering strengths much smaller than the subband width, 
and an asymmetry in the component density of states. 

The most notable result is, however, the large en- 
hancement of the spectrally resolved return probability 
P(tt', 77) for small but finite 77 in the narrow energy re- 
gions where the dispersion of the polaron states becomes 
extremely flat. This occurs in the absence (presence) of 
alloy-type disorder at the high-energy edge(s) of the po- 
laronic subband (mini-subbands) . To elucidate the phys- 
ical meaning of this strong enhancement, we analyzed 
the 77— asymptotics of P{Lu,ri) and introduced the con- 
cept of a spectrally resolved delocalization time which is 
the characteristic time scale on which, at a given energy, 
the electron leaves a given site. The strong enhance- 
ment of P{uj, rj) in the flat part(s) of the subband (mini- 
subbands) signals large delocalization times. According 
to their delocalization times we could thus classify, within 
the limitations of the DCPA, "fast" quasi-particle-like 
polaron states at the bottom and "sluggish" defect-like 
polaron states at the top of the subband (mini-subbands). 

Within the DCPA, we could not decide, however, 
whether "defect-like" , i.e., temporarily trapped states are 
in fact "defect", i.e., trapped states. At this point it is 
important to recall that within the DCPA only on-site 
processes contribute to the delocalization times. In con- 
trast, inter-site processes are completely neglected. Par- 
ticularly, back-scattering processes are expected to have a 
significant effect. In the conventional alloy problem, for 
example, these processes even yield, depending on the 
model parameters, infinite delocalization times at cer- 
tain energies. It is therefore conceivable that the DCPA 
results for the delocalization times of polaron states pre- 
sented in this paper are only lower bounds and temporar- 
ily trapped states may in fact be trapped states. More 
sophisticated methods beyond the DCPA, taking back- 
scattering and other non-mean field effects explicitly into 
account, are clearly needed to further elucidate the trap- 
ping issue. 
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APPENDIX A: POLARON IMPURITY MODEL 
(PIM) 

In this appendix we show how the local two-point func- 
tion Gii{z) can be calculated from the PIM. Towards that 
end, we apply < i\{N,ni\...\mi, N)\i > on both sides of 
Eq. (p7|) and obtain in compact matrix notation 



I?W(z) = GW(z) + GW(z) 

'ph 



with matrices 

2?SW = <*l(0,n,|i?f)(z)|m„0)|*>, 



(Al) 

(A2) 
(A3) 



gVn+ l(5,„,„+i - gVnSm.n-i, (A4) 

(A5) 



To derive Eq. ( |Al| ) we employed AiJi(z) ^ \i >< i\ to- 
gether with G^'^\z) —< i\g'^'^\z)\i >, which follows from 
the definition of Geff{z). 

If we now introduce an auxiliary matrix 

hW(z)= [i-GW(z).W(5)1"'gW(z), (A6) 



which sums all electron-phonon scattering processes, Eq. 
(Al) can be rearranged into 



pW(z)=[l-ifW(z)ai:V..;z)l" i/W(z). (A7) 



The configuration average is straightforwardly per- 
formed since (single-site) randomness enters only through 
o'i^pixi] z). We find, using the bi- modal probability dis- 
tribution p{xi) — cS{xi) + (I — c)S{xi — 1), 



((I?(~)(z))) - J dx,p{x,)V^''\z) 

= il~c) [l-i/(^n-)Ol;-) 
X #W(z) 

+ c[l-#W(z)ai^^(0;z)]'' 
X ijW(z) 

= (l-c)I?^^W(z) + cI?^^W(z), 
or, explicitly in terms of matrix elements, 
((pW(z))) = (1 - c)I?;^;W(z) + cl?fi^)(z). 



(A8) 



(A9) 



where the amplitudes "Dni^^ {z) satisfy a set of recursion 
relations (A = A, B) 



V^^^\z) = F^'(^+")(z)5„,™ - gF^-^^^+-\z) 
with 

F^'(")(z) = 

1 

z — — ex + Gii^{z - nO) - TZ[Qii{z — nfi)] 



(AlO) 



(All) 



From the self-consistency condition, Eq. (29), finally 
follows 



fc(z) = ((I?^"o)(z))) 

= (l-c)I?^^(°)(z) + cPo^„'(")(z), 



(A12) 



which directly yields Eq. ( pO| ) if the continued fraction 
expansion |Q for the amp litudes 'Dqq^^ (z) is invoked to 



iteratively solve Eq. ([A1C|) 



APPENDIX B: DERIVATION OF EQ. (g) 

In this appendix, we bring Eq. (p3|) into the numeri- 
cally more convenient form of Eq. (p4|). As a preparatory 
step we notice, using the definition of the PIM, Eq. ([2^), 
that the matrix elements of the atomic T-matrix defined 
in Eq. (56|) can be written as 



Gtiiz)gii{z - qQ) 



(Bl) 



with V^J{z) defined in Eq. (A2); here, we employed 

again the identity Gii \z) —< i\g^'''' {z)\i >. If we now 
insert Eq. (Bl) into Eq. ( p3| ) and take into account that 

((I?(°)(z))) = Quiz), Eq. iM) reduces to 



OO 



/(zi,Z2) = }_^{{V^^Jiz^)V'^'iz2))). (B2) 

q=0 



For zi 



z*, we may use z)'-*'^(z) = [7j("^(z*)]* and 



(0)( 



slightly rearrange the above equation into 



fiz*,z) 



q=0 



(B3) 



Employing now the matrix notation introduced in ap- 
pendix A, the configuration average of Eq. (B3) can be 
performed explicitly. To be specific, from Eq. (A7) it 
follows that 



v: 



X Hl^\z) 



(B4) 



which, using the bi-modal probability distribution 
p{xi) = c6(xi) -t- (1 — c)6(xi — 1), immediately yields 
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(l-c)II?, 



qO 

qO 



c\V. 



B,iO) 
qO 



(B5) 



with the amphtudes Ty^Q ^"^ (z) satisfymg the recursion re- 



lations Eq. ( |A1C| ). Combining Eqs. ( |B3| ) and ( |B5[ ) and 

setting z — Lu + ir], we finally obtain Eq. (54). 
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FIG. 1. Schematic representation of a fourth order pro- 
cess contributing to Ql . Atomic T-matrices and effective 
one-resolvents are, respectively, depicted by solid circles and 
solid lines. The arrows indicate the order in which the dia- 
gram has to be read. The energy arguments are clear from 
the context. 



(a) 



(b) 



.,C---Q '" 



FIG. 2. Schematic representation of the second order (lad- 
der) process contributing to the local DCPA vertex (a) and 
of the second order (maximally crossed) process neglected 
within the DCPA (b). The dashed line depicts the operator 
Oe and solid circles and lines stand, respectively, for atomic 
T-matrices and effective one-resolvents. The arrows indicate 
the order in which the diagrams have to be read. The energy 
arguments are clear from the context. 



+ 



+ 



FIG. 3. Schematic representation of Eq. (^0|). The box 
on the left hand side denotes the local vertex operator Ti, 
the dashed line depicts the operator Oe, and solid circles and 
lines stand, respectively, for atomic T-matrices and effective 
one-resolvents. The arrows indicate the order in which the 
diagram has to be read. The energy arguments are clear from 
the context. 
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FIG. 4. Spectral and diffusion properties of the ordered 
Holstein model in the vicinity of the first polaronic subband 
for A = 1.96, a = 2.8 [see Table |, and = 10"". (a) Density 
of states N{uj) (sohd line) and imaginary part of the coher- 
ent potential Irav(u) (dashed line), (b) Spectral function 
A{ijj,e^) with eg ranging from eg = —1 to eg = 1 (bottom 
to top) in steps of Aeg = 0.2. For clarity the spectral func- 
tions corresponding to different eg are artificially shifted along 
the vertical axis, (c) Spectrally resolved return probability 
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FIG. 5. Blow-up of Fig. ^ in the immediate vicinity of the 
high-energy edge of the lowest polaronic subband. (a) Real 
part of the coherent potential Rev{Lo) for rj = 10~* (thick solid 
line). The intersections (below lo — —1.125) with the thin 
solid and dashed lines determine, respectively, the position 
of the maximum of the subband density of states (wm) and 
the position of the high-energy edge (tOh)- (b) Density of 
states for 77 = 10~^ (solid line) and rj — 10~* (dotted line). 
The dashed line depicts the imaginary part of the coherent 
potential Imv{uj) for rj — 10~^. (c) The poles of the spectral 
function io{ej:). Here we plot ej: vs. uj. (d) Spectrally resolved 
return probability PioJ^rf) for 77 = 10"*. 
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FIG. 6. Spectral and diffusion properties of the Holstein al- 
loy in the vicinity of the high-energy edge of the first polaronic 
subband for c = 0.5 and 5 = —0.01. The polaron parameters 
are A = 1.96 and a = 2.8 [see Table | and r; = 10"*. (a) Den- 
sity of states N{uj) (solid line), component density of states 
Na{u}) (dotted line) and Nb{u}) (dot-dashed line), and imagi- 
nary part of the coherent potential Imv{uj) (dashed line), (b) 
Spectral function y4(tj,eg) with e^: ranging from eg = —1 to 
eg = 1 (bottom to top) in steps of Aeg = 0.2. For clarity the 
spectral functions corresponding to different eg are artificially 
shifted along the vertical axis, (c) Spectrally resolved return 
probability P{uJ,ri). 
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FIG. 7. Density of states N{uj) for the Holstein alloy in the 
vicinity of the first polaronic subband. The polaron param- 
eters are A = 1.96 and a — 2.8 [see Table Q and rj — 10~*. 
The scattering strength 5 — —0.004 and the concentration of 
the B atoms varies from c — 0.0 (bottom) to c = 1.0 (top) 
in steps of Ac = 0.2. For clarity the density of states cor- 
responding to different concentrations are artificially shifted 
along the vertical axis. 
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FIG. 9. ry— asymptotics of the spectrally resolved return 
probability P{ij,ri) for the ordered Holstein model. The 
polaron parameters are A — 1.96 and a — 2.8 [see Table 
|]. The three sets of curves correspond to three energies: 
LO = -1.135,-1.127, and -1.12549 (from bottom to top). 
Solid [dashes] lines depict the numerical data [asymptotic ex- 
pansion given in Eq. (|60|)]. The corresponding density of 
states for 77 = 10~* is shown in the inset. 



CL 0.5 




CO 

FIG. 8. The lower panel shows the density of states N{uj) 
and the imaginary part of the coherent potential iTnv(u) for 
the Holstein alloy in the vicinity of the high-energy edge of 
the first polaronic subband for c = 0.5 and three different 
values of the scattering strength 5. The polaron parameters 
are A = 1.96 and a = 2.8 [see Table § and 77 = lO""*. The 
corresponding spectrally resolved return probabilities P{w, rj) 
are given in the upper panel. 
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FIG. 10. Spectrally resolved delocalization time To{uj) for 
the ordered Holstein model in the vicinity of the high energy 
edge of the first polaronic subband. The polaron parameters 
are A = 1.96 and a = 2.8 [see Table ^. The inset shows the 
density of states for rj = 10~*. Note, the data for To(tj) are 
for energies below the rj ^ high-energy edge. 
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FIG. 11. Spectrally resolved delocalization time To{uj) for 
the Holstein alloy in the vicinity of the high energy edge of 
the first polaronic subband for c — 0.5 and 5 = —0.002. The 
polaron parameters are A — 1.96 and a = 2.8 [see Table |. 
The inset shows the density of states for rj = 10~*. Note, the 
data for ro(tj) are for energies below the r; — > high-energy 
edge. 



TABLE I. Polaron model parameters (in units of 2J) 



J 


n 


g 


A = 2g'^/n 


a = g/n 


0.5 


0.125 


0.35 


1.96 


2.8 
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